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This paper investigates the effect of nonlinear large deflection bending on the aerody- 
namic performance of a high aspect ratio flexible wing. A set of nonlinear static aeroelastic 
equations are derived for the large bending deflection of a high aspect ratio wing structure. 
An analysis is conducted to compare the nonlinear bending theory with the linear bending 
theory. The results show that the nonlinear bending theory is length-preserving whereas 
the linear bending theory causes a non-physical effect of lengthening the wing structure 
under the no axial load condition. A modified lifting line theory is developed to compute 
the lift and drag coefficients of a wing structure undergoing a large bending deflection. 
The lift and and drag coefficients are more accurately estimated by the nonlinear bending 
theory due to its length-preserving property. The nonlinear bending theory yields lower 
lift and span efficiency than the linear bending theory. A coupled aerodynamic-nonlinear 
finite element model is developed to implement the nonlinear bending theory for a Com- 
mon Research Model (CRM) flexible wing wind tunnel model to be tested in the University 
of Washington Aeronautical Laboratory (UWAL). The structural stiffness of the model is 
designed to give about 10% wing tip deflection which is large enough that could cause the 
nonlinear deflection effect to become significant. The computational results show that the 
nonlinear bending theory yields slightly less lift than the linear bending theory for this 
wind tunnel model. As a result, the linear bending theory is deemed adequate for the 
CRM wind tunnel model. 


I. Introduction 


Modern aircraft are increasingly designed to be highly aerodynamically efficient in order to reduce the 
fuel consumption, hence operating cost of modern transports. To achieve this goal, the aircraft industry has 
been adopting a high aspect ratio wing design with the advanced composite construction in modern trans- 
ports. Aircraft design concepts that take advantage of wing flexibility to increase aerodynamic performance 
and maneuverability have been investigated. By twisting a wing structure, an aerodynamic moment can be 
generated to enable an aircraft to execute a maneuver in lieu of the traditional flight control surfaces. For 
example, a rolling moment can be induced by twisting the left and right wings in the opposite direction. 
Similarly, a pitching moment can be generated by twisting both wings in the same direction. Wing twisting 
or warping for flight control is not a new concept and was used in the Wright Flyer in the 1903. The U.S. Air 
Force conducted the Active Flexible Wing program in the 1980’s and 1990’s to explore the potential use of 
leading edge slats and trailing edge flaps to increase the control effectiveness of F-16 aircraft for high speed 
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maneuvers.! In the 2000’s, the Active Aeroelastic Wing research program also investigated a similar wing 
twisting technology for improved roll maneuverability of a modified F/A-18 aircraft.2, Wing shaping control 
concepts for drag reduction are being studied by NASA to leverage wing flexibility for aerodynamic per- 
formance.?}* By re-twisting a flexible wing and using variable camber aerodynamic flight control surfaces, 
aircraft wings can have a mission-adaptive capability to optimize L/D throughout a flight envelope.? In 
recognition of the role of aeroelasticity on aircraft performance and dynamics, NASA Advanced Air Trans- 
port Technology (AATT) project is conducting research in the area of Performance Adaptive Aeroelastic 
Wing (PAAW). This research develops concepts such as the variable camber continuous trailing edge flap 
(VCCTEF) to enable wing shaping control for aerodynamic performance and dynamics.**® Currently, a 
wind tunnel test is being planned to conduct an active real-time drag minimization control strategy for 
a high aspect ratio flexible wing based on the Common Research Model” equipped with the VCCTEF in 
the University of Washington Aeronautical Laboratory (UWAL). The wing is purposely designed with a 
reduced structural stiffness to achieve about a 10% wing tip deflection. Efforts are under way to create an 
aerodynamic-structural model of this wind tunnel model. It is of interest to investigate the effect of nonlinear 
large deflection bending on this wind tunnel model so that a better understanding can be ascertained prior 
to the wind tunnel test entry. 

Structural deflections of lifting surfaces interact with aerodynamic forces to create aeroelastic coupling 
that can affect aircraft performance. Understanding these effects can improve the prediction of aircraft 
performance and provide an insight into how to design an aerodynamically efficient high aspect ratio flexible 
wing. Generally, high aspect ratio lifting surfaces undergo a greater degree of structural deflections than low 
aspect ratio lifting surfaces. As a result, the nonlinear structural deflection could have a significant effect on 
the aerodynamic performance. 

The purpose of this study is to investigate the effect of nonlinear large deflection bending on the high 
aspect ratio CRM flexible wing to ascertain whether or not this effect should be given further consideration 
in the analysis of the wind tunnel test data and future optimization efforts. Toward this end, a nonlinear 
bending theory is developed along with a modified lifting line theory which includes the effect of the wing 
aeroelastic deflections. A nonlinear finite element model is developed to implement the nonlinear bending 
theory. The finite element model is then coupled to an aerodynamic model of the CRM wing which is based 
on a vortex lattice method with transonic and viscous corrections via a 2D transonic small disturbance 
method? and a 2D integral boundary layer method.!° The coupled aerodynamic-nonlinear finite element 
model is analyzed to quantify the aerodynamic performance of this model.® 


II. Common Research Model 


The aerodynamic model is a wind tunnel sub-scale model of a CRM wing to be tested in the UWAL in 
2017 to demonstrate an active real-time drag optimization control strategy. This CRM model is designed to 
have about 10% wing tip deflection to represent the current state of the art high aspect ratio wings in modern 
transport aircraft such as the Boeing 777. The model is 85 inches in length with an elastic axis sweep angle 
of 31.5° and includes six active two-segment VCCTEF control surfaces as shown in Fig. 1. These surfaces 
are driven by actuators for the real-time drag minimization control strategy. The test conditions are at low 
subsonic speed of about Mach 0.1. 

An aerodynamic model of this CRM wing is constructed using NASA vortex lattice code VORLAX. 
A transonic correction has been developed for VORLAX to account for transonic flow over wings. The 
transonic correction is based on the transonic small disturbance (TSD) theory® and is implemented using a 
2D transonic code TSFOIL developed by NASA Ames Research Center.? For the wind tunnel model analysis, 
the transonic correction is not needed, but for a full-scale CRM model, this correction will be required. In 
addition, the viscous effect is also corrected for in VORLAX via a 2D integral boundary layer code developed 
in-house by NASA Ames Research Center.1? The coupled transonic plus integral boundary layer method 
with the vortex lattice method has been shown to be able to produce reasonably accurate aerodynamic 
performance prediction as compared to a Reynolds-Averaged Navier-Stokes (RANS) CFD solver which is 
computationally much more expensive than the vortex lattice method.? Moreover, the vortex lattice code 
VORLAX has the mesh deformation capability that allows it to be easily coupled to a finite element model for 
static aeroelastic analysis,'! whereas the mesh deformation capability is only available for a few NASA RANS 
CFD solvers. With the nonlinear finite element capability, the coupling will require more iterations for a 
solution to converge. Thus, this will create an additional computational cost that must be taken into account 
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for RANS solvers, whereas such cost is very small for the vortex lattice method even with the transonic and 
viscous corrections. For conceptual aerodynamic analyses and multidisciplinary design optimization, this 
method is very computationally efficient and can generate rapid candidate design solutions. 


Figure 1. Common Research Model Sub-Scale Wind Tunnel Model in UWAL 


Figure 2 shows the surface pressure distribution of the full-scale CRM wing at Mach 0.85 and Reynolds 
number of 1.88 million per foot. The model shows a transonic region in the wing root area as is typically 
the case with transport wings. Figures 3 and 4 shows the computed lift curve and drag polar of the CRM 
rigid wing. 
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Figure 2. Surface Pressure Distribution of CRM Rigid Wing Model at Mach 0.85 Computed by VORLAX 
with Transonic and Boundary Layer Corrections 
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CRM Rigid, Mach=0.85, Re=1.88x10°/ft 
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Figure 3. Lift Curve of CRM Rigid Wing Computed by VORLAX with Transonic and Boundary Layer 
Corrections 
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Figure 4. Drag Polar of CRM Rigid Wing Computed by VORLAX with Transonic and Boundary Layer 
Corrections 


The aeroelastic effect on the aerodynamic performance is typically described by the following equation 
for the local aeroelastic angle of attack defined with respect to the elastic axis according to the linear bending 


theory: 
a 
Oe = oq We tanA—© (1) 
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where & is the coordinate of the elastic axis, a is the wing angle of attack, A is the sweep angle of the elastic 


axis, © is the torsional twist (positive nose-down), W is the bending deflection, and the subscript x denotes 
(2) 


. Thus, the wing bending deflection generally causes a reduction in lift due to the decrease in the local 
aeroelastic angle of attack for a sweptback wing with A > 0. The torsional twist is generally nose-down 
(positive in the sign convention). Thus, the net effect is a reduction in lift as well as a change in the design 
spanwise lift distribution which can affect the induced drag. If the effect of the nonlinear large deflection 
bending is significant, then clearly this cannot be neglected in the aerodynamic performance analysis of a 
high aspect ratio flexible wing structure. 


III. Nonlinear Large Bending Deflection Analysis 


To analyze the nonlinear large deflection bending effect, we perform a strain analysis. Toward that end, 
consider an airfoil section on the left wing as shown in Fig. 5 undergoing a combined axial, bending, and 
torsional displacement field. Let (x,y,z) be the undeformed coordinates of point Q on a wing airfoil section 
in the reference frame D defined by unit vectors (di,d2,d3). Let po = xdj be a position vector along the 
elastic axis. Then, point Q is defined by a position vector p = po + q where q = ydz + zd3 defines point Q 
in the y — z plane from the elastic axis. Then, the undeformed local airfoil coordinates of point Q are 


y |__| cosy —siny ff] (2) 
z siny cosy é 
where 7 and € are local airfoil coordinates, and ¥ is the wing section pre-twist angle, positive nose-down. 
Differentiating y and z with respect to x gives 


foe oir] L a | : 
2a cosy —siny € yy 
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Figure 5. Left Wing Reference Frame of Wing in Combined Bending-Torsion 


Let O be a torsional twist angle about the z-axis, positive nose-down. Let W and V be flapwise and 
chordwise bending deflections of point Q, respectively, and W is allowed to be large relative to V. Let U be 
the axial displacement of point Q. Then, the displacement and rotation vectors due to the elastic deformation 


can be expressed as 
r = Ud, + Vd2+ Wd3 (4) 


o = Od, — sin™' W,d2 + V,d3 (5) 


where the subscripts x and s denote the partial derivatives of W and V with respect to x and the arc length s. 
Note that the term sin~' W, represents the rotation angle due to the large flapwise bending with respect to 
the arc length. If W, is small, then the usual approximation sin~' W, + W, applies. Otherwise, W, = W.sz 


applies where 
ds \2 
se = = VIF BH Vt +2) 7) (6) 
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Let (21,41, 21) be the deformed coordinates of point Q on the airfoil in the reference frame D and 
Pi = 21d; + y1dg + 21d3 be its position vector. Then the coordinates (21, y1, 21) are computed as 


pi=p+r+¢xq (7) 
where 
Ly z+U —yV, — zsin ! W, 
yo | = yt+V—-2z0 (8) 
zy z+W+y0 


Differentiating x1, y1, and z, with respect to x yields 


1+ Uy, — yVrx 4 zy Ve z—Wesse yy sin” ' W, 


U1 1-W2 
Ya | = zy + V, — 20, —yy 8 (9) 
Z12 yy +sin-!W, + yO, — 27 ® 


Note that the term eae represents the physical curvature of the nonlinear bending with a large 


deflection. This term is consistent with the nonlinear bending curvature formula by Hodges.!* Other 


authors!®:!© have used the curvature formula Fi ee iH based on the elementary calculus which is sometimes 
+W2)2 
mistaken as the physical curvature of a beam as stated by Hodges.'4 The curvature formula tes can be 


s 


used in the undeformed coordinates whereas the curvature formula ae must be used in the deformed 
14+W2)2 


coordinates which are impractical for implementation. 
Neglecting the transverse shear effect, the longitudinal strain is computed as! 
ds,;—ds_ s 
_dey—ds _ sie | (10) 


ds Sx 


where 


— 2 2 2 
Sia = Vite te Ye + cae 


= 2 + 2U, — 2yVrw — 22 


$1,” is approximated by a Taylor series expansion as 


2 1\2 1\2 

Wate cag arate, « Ce Dt (metet)' + (He) 
+(y? +27) 7,4 tee (12 
coy tw te) (12) 


s 


Sia LY Sz Uz YVrw 


The slope of the twist angle 4’ can play a significant role in structures with large twists such as turboma- 
chinery blades. For aircraft wings, this effect is negligible and therefore can be neglected. Thus, for a small 
wing twist angle slope, + ~O0and s, +1. Then W, % W, and W,, © Wee. Thus, the longitudinal strain is 
then obtained as 


_ Were PT ae ee Le ee ee 
€ = Uz — yen — 2 ; Wat v2 + Ve t 5 (sin W,) +5 \(y + 2?) 03 + 59°Vee 
1 2 Wre Were \ Wee sci 
ae T— we | YUeVen 2Uy Tw + yz Ving Tw zV,0, +ysin-°W,0, (13) 


Assuming that the neutral axis of a wing section passes through the elastic axis, then the axial force and 
moments acting on a wing are evaluated as® 


1 1 1 
P= , EedA = EAU, + 5EA lv? + V? + (sin7* W2)’| + 5 Flex; + 5 EliaVee 
il w2 W. 
= El] yy,—*2. — | EAegUy + Elyz,———“— | Vox (14 
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M, = GJO, + i Ee (y? + 2°) Q,dA = 


Were 
CIP EL SaeD (20. a BB. | 


1-W2 
1 Brat 1 w2 W. 
~E lye [v2 y2 : -ly, - EBwW?2 “BB. **_ + FEB Vo = - 
ae 24+V2 + (sin ) + 5 EBWir + 5 5y_wat 6 Ww 
1 
+ (-EB2V, + EBssin“'W.) ©; + 5EBi@; (15) 
Wee 1 ae : 
M,=— | EezdA=(1+U;) { Elyy =a ElyzVex | — 5BB20, — 5 EBV er 


W. Woe a 
EBs —+, — EBoVix <=> + (ElyyVz + Elyzsin™'W,) Ox (16) 


Ji-W3 


Were 
1-W?2 


x 


1 1 
+ E..¥e) 7 5 EB38, - 9 EBV ee 


W. Wrz Sees 
EBg ae EBV." € (ElyzVe + EI,, sin : We) Ox (17) 


Vi-Wi 


where E£ is the Young’s modulus, G is the shear modulus, A is the tensile area, J is the torsional constant, 
and B;,i=1,2,...,10 are defined as 


By, (y? + z7) 

Bo Zz (y? + a) 

Bs y(yrt2? 

B, y? (y? L 2) 

Be ||: 2 (y et: 2) 

re | = / vette |e (18) 
Br Yyrz 

Bg 3 

Bog ye? 

Bio y° 


For a rectangular wing box structure with the neutral axis in coincidence with the elastic axis, the 
constants B,, B,, and Bs are generally much larger than the other constants B;, i 4 1,4,5 and can be 
approximated as 


B mee 
By ~ A Lealzz (19) 
Bs Lealyy 


Assuming that the nonlinear contributions of the chordwise bending deflection are small and neglecting 
the cross-product inertia J, and the small constants B;, i 4 1,4, 5, then the following simplification can be 
made: 


P, = EAU, + +BA [v2 + (sin! W, )’) + 1 Blge92 + SEL We (20) 
M, = (GJ + Port) Oz (21) 

Wes 
M, = ElzVex — Elz, sin” | Wz9z (23) 


where rz = 4/ fe is the radius of gyration. 
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Note that the signs of the moments are defined in the 


positive deflection sense such that 


M = M,d, — M,d2 + M,d3 (24) 
The resulting equilibrium equations are given by® 
or Py ty (25) 
Pe ce 
9 M, PyW, — PzVz — Mz 
ag | Me | =| PeWe-P.+my (26) 
M, P,Vx — Py -— mz 
These equations become 
OM, OM, OM, 
+ P,Vy +m.) We “— P,W, —my ) Vi +m, =0 (27) 
Ox Ox Ox 
0 0 
O?M, O(P:Wz) Omy 
20 2 
Ox? Ox Ox f2=0 2) 
O°M, O(P,V,) . Om, 
=0 29 
Ox? Ox Ox ty ey) 
The resulting nonlinear equations are obtained as 
© {pau,+iBA [U2 + (sin? W)"] Aes cup Were ong (30) 
Ox NG 7 is QE OE We 7 
a 2 
o? Wea O(P,W,) om 
—, |14+U,) EL =f. y 32 
Ox? (1+ Uz) Fly 1-W2 Ox fa+ Ox (32) 
OP O(P, Vx) Om, 
as (ElzzVex — Elz, sin~' W292 = 33 
Ox? ( ia ) Ox ty Ox of 


IV. 


The forces and moments on the right hand side of Eqs. 


Static Aeroelasticity 


(30) - (33) include all external and inertial forces 


and moments due to the acceleration of the wing structure. In order to compute the aerodynamic forces and 


moments of a flexible wing structure, static aeroelasticity 


must be considered. 


The relative velocity of the air approaching a wing section includes the contributions from the wing 
aeroelastic deflections that result in changes in the local angle of attack. Since aerodynamic forces and 
moments are dependent on the local angle of attack, the wing aeroelastic deflections will generate additional 
elastic forces and moments. The local angle of attack depends on the relative approaching air velocity as well 
as the rotation angle @ from Eq. (5). The relative air velocity in turn also depends on the deflection-induced 
velocity. The velocity at point Q due to the aircraft velocity and angular velocity in the reference frame D 


is computed as 


vg =V+w xr = (ub) + vb2 + wh) + (pbr 


qb2 + 1rb3) X (—2abi — Yabo — Zab3) 


= (u +TYa — qa) by + (v — Tha + pza) be 4 
= ©, + yrd2 + %d3 
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where (b;, bz, bs) are the unit vectors in the standard aircraft body-fixed reference frame B. The coordinate 
transformation between the reference frame B and D is given by 


b; —sinAcosI’ —cosA sinAsinI dy 
bz | =| —cosAcos[f sinA  cosAsinI dy (35) 
bs —sinT 0 —cosT ds 


where I is the wing dihedral angle. 
(2, Yt, 24) are the velocity components in the left wing coordinate reference frame D and are computed 


as 
Lt —(utrya — Ga) sin AcosT — (v — rag + pzq) cos AcosT — (w + q&a — pya) sinT 
4 |= —(ut+rya — qa) cos A+ (vu — rag + p2q) sin A (36) 
Zt (u+rya — Q2a) sin AsinT + (v — rag + pz) cos AsinT — (w + qa — pya) cosT 


where (u,v, w) are the aircraft velocity components, (p,q,1r) are aircraft angular velocity components in the 
roll, pitch, and yaw axes, and (2a, Ya, Za) is the coordinate of point Q in the aircraft body-fixed reference 
frame B relative to the aircraft CG (center of gravity) such that x2, is positive when point Q is aft of the 
aircraft CG, ya is positive when point Q is toward the left wing from the aircraft CG, and z,q is positive 
when point Q is above the aircraft CG. 

Consider a trim problem when 6 = 0, p = q = r = 0 for steady-state aerodynamics. For simplicity, we 
assume I = 0. Then, 


Lt —usin A 
ye | = | —ucosA (37) 
Zt —W 


The velocity must be transformed from the reference frame D to the airfoil local coordinate reference 
frame defined by (1,7, €) as follows: 


Un 1 0 0 cosV,  sinV, 0O 
Un | = | 0 cos sinO —sinV, cosV, O | X 
VE 0 —sin© cosO 0 0 1 
cos (sin~'W,,) 0 sin (sin W,) Ly 
x 0 1 0 Yt (38) 


— sin (sin7! W;) 0 cos (sin7* Wz) 2t 


Making the small angle assumption for O and V, only and neglecting the nonlinear terms that contain 
V,, then the velocity components in the local coordinate reference frame are computed as 


Dis tef1—W2 + yVe + We 


vy | =| —2t (Ve + We) + ye + 4 OV/1— W2 (39) 
VE tiW yO b+ Z¢4/ 1 Ww2 


For a small angle of attack, the usual assumption of a = @ is valid, but if the rotation angle due to 


the flapwise bending is large, then some accuracy in the small angle of attack assumption may be suffered. 
Therefore, the exact angle of attack definition a = tan~ * is used. Then, the local aeroelastic angle of 
attack on the airfoil section due to the velocity components v,, and ve defined with respect to the elastic axis 


is computed as 


v ks 
a, = tan? “5 w tan“! 6 4 "5 (40) 
aaarTc) 
where 
V¢e =—-W (41) 
Ave = usin AW, + ucos AO — w (/1- W? - 1) (42) 
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Uy, = —ucos A (43) 


Av, = usin A (V, + OW,) — wOvV1—- W? (44) 
Then, 
" 1 wsin AW, + weos AO — w (/1 = W? - 1) 
Q, = tan™+ ( ) 7 
ucos A 1+ 545 ucos A 
1 —uw sin A (V, + OW,) + w?9,/1 — W2 (45) 
l+ oo we u? cos? A 


2 
Using the Taylor series approximation of ,/1 — W? + 1— We and upon simplification, we get 


2 
tan-2 ( 222 cos? A sin AW, + cos AQ + oWs 
= ta 
re ” cos A cos? A + tan? a cos A 
oe 
cos? A —tanasinA(V, + OW,) + tan? aQ (1 _ We ) 
cos? A + tan? a cos? A (46) 


Using the Taylor series approximation of tan7' (#29) ~ ‘ne (1 jut a.) mw (1+ «) 5 


cos A cos A 6 cos? A cos A 


cos A 


2 
. f oes x (1 =) ~~ (1 - on tan? A), the local aeroelastic angle of attack can be expressed as 


Ac = Arp + Ae (47) 
where Bea 
a a’ tan 
= 1 48 
Or cos A ( 6 ) es) 
is the rigid angle of attack due to the aircraft velocity and 
tan asin AV, cos A — tanaO tanaW2 
— in AW, + ——— ] -09 49 
ve cost A+ tan2a cos? A+ tan2a (sin zs ) a) 


is the elastic angle of attack due to the aeroelastic deflections. Note that the rigid angle of attack is defined 
with respect to the elastic axis, hence the term cosA in the denominator. For a sweptback wing, the 
contribution of W, is negative, thus effectively reduces the local angle of attack. The bending deflection thus 
creates an effective wash-out twist to reduce the local angle of attack for sweptback wings. Consequently, 
the trim angle of attack must be increased in order to compensate for the wash-out twist. Equation (49) for 
nonlinear large deflection bending should be contrasted with Eq. (1) for linear small deflection bending. 


A. Forces and Moments 


The aerodynamic and body forces and moments are given by 


te Cpqdooc sin A cos A 

i CDdooC cos? A 

fe = CL dooc cos A — mg (50) 
My —CmooC? cos” A + Mgecq 

My Cm oo sin A cos A 

Mz 0 


where qo is the dynamic pressure, c is the chord length, m is the mass distribution, g is the gravity 
acceleration, and €¢, is the offset of the center of mass from the elastic axis (positive with the center of mass 
aft of the elastic axis). 
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The steady-state lift, drag, and pitching moment coefficients in the streamwise direction can be generally 
modeled as 


CL = Cro + Cx, (@) a. cosA (51) 
cp = Cp, + ep, (a) a cos A + ep2 (a) a% cos” A (52) 
Ci Cae (Gl beh (53) 

Cc 


where cy, and cp, are the lift and drag coefficients at zero angle of attack, cy, is the lift curve slope, cp, 
and cp_, are the drag derivatives with respect to the angle of attack, c,,,, is the pitching moment coefficient 
at the aerodynamic center, and e is the offset of the aerodynamic center from the elastic axis, positive with 
the aerodynamic center forward of the elastic axis. The nonlinear effect of large angles of attack in the 
aerodynamic coefficients is reflected in the general functions cz, (a), cp, (@), Ep, (@), and Cm,, (a). These 
functions take on constant values for small angles of attack below stall. 


B. Effect of Nonlinear Bending Theory on Lifting Line Aerodynamics 


Consider the case of unswept wings, then f, = 0. The local angle of attack is reduced to 


1—tana® tanaW?2 
1+ tan? a 2 


Ae =a (S) (54) 

For linear structures, the wing bending deflection has no effect on the wash-out twist for unswept wings. 
In contrast, for nonlinear structures, the contribution of the wing bending deflection to the wash-out twist 
is non-zero due to the term W2, although this term is small. Nonetheless, this results in a bending-torsion 
aeroelastic coupling even though the wings are unswept. This coupling does not exist in linear bending 
theory. This term can significantly affect the lift distribution and the induced drag of a high aspect ratio 
flexible wing as will be shown. 

Since there is no applied axial force at the end of the cantilever wing, the axial displacement equation is 
equal to zero 


P, = EAU, + 3EA [U2 + (sin-! We)"] + 1 i tpeO? + LBIy 22, =0 (55) 
Solving for U, from Eq. (55) yields 
_ omnes | 2 Low Lyy W2, 
1+U, = [t= Gon W;) A 2 AT-W2 (56) 
Then, the nonlinear static aeroelastic equations become 
3) 1—tana® tanaW? é 
Du (GIO@z) +er. Gee ( (Stata 5 °) = Cm,.qooC” — MgCcg (57) 
in 2 Lyy W?, 
2 1 — (sin“' W,)" — r202 — 4 22, 1 —tana® tanaW?2 
as | ElyyWee = 1 +¢7 doo€ 7 +0 
Ox? 1-W? ° 1+tan*a 2 
=Cr,.dooc— mg (58) 
oP “ease 1— tana® tanaW? 
B22 (Elz2Ven — El,, sin W,9-2) + (co, + 2cp2 a) doo€ ( aes 5 + °) 


1—tanaO tanaW2 : 
oa Be Z—O] = Pe 59 
ears ( 1+tan?a 2 ) CDrdooe (59) 


where cp, = cr, + cx, is the rigid lift coefficient, cp, = cp, +cp,a+ CD.2 a? is the rigid drag coefficient, 
and Cm, =Cm,. + x, © is the rigid pitching moment coefficient about the elastic axis. 
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If the torsional twist is small relative to the flapwise bending slope, then the nonlinear static aeroelastic 
equation for bending is expressed as 


Cas 2 Iyy W2, 
ci ElyyW. pe i ee IeWa | aie , ee =C c—m (60) 
Ox? mead 1-W2 Boa doo%5 (1+tan?a) | ls : 


For the nonlinear large deflection bending of the wing structure, a modified lifting line theory is used. 
The velocity induced by a trailing vortex at x = xo acting at the quarter chord is given by the Biot-Savart 
law according to!” 


rx V ad 
War AnVoV.r ds (61) 
where ae is the lift circulation distribution along the arc length s of the wing, r = — (a — «+ Up — U) d, — 
(Wo — eOo — W + €O) dz, and V = Vida. This yields the following expression: 
d - - d, dl 
dv; = (xo x Uo U) 3+ (Wo €Oo W+ eO) 1 d as (62) 


Arrv.r ds 


Due to the axial displacement and flapwise bending effects, in addition to the induced downwash velocity 
component, there also exists an induced sidewash velocity component. The effective downwash along the arc 
length of the wing is the induced velocity component normal to the vortex sheet.!’ This is given by 


dw; = dv;.No (63) 


dso 


where No = [W204 = (1 + Uz0) d3| das and dso = dro y/(1 + Ueal + W2 0. Then, 


d= (to — & + Up — U) (1+ Uso) + (Wo — G0 — W + €©) Wao dxo dD (64) 


An (xo —~2£+Up—U)? + (Wo — e@y — W + €0)” dso ds 


where 09, Uo, Uz, Wo, and W;z,o are the wing torsional twist, axial displacement, axial displacement 
derivative, bending deflection, and bending slope, respectively, at x = 20. 
The total induced downwash at x = xo is then evaluated by the following expression: 


b 
2 - Uj —-U)(1+U, Wo — eO0 — W + €9) Wz 9 dx aT 
wi (xo) = i, vores LAL o) + ites aie ) On de (65) 
-3 An (xo x+ Uo U) t (Wo eOo W 4 eO) dso dx 
Note that if U = 0 and U, = 0, Eq. (65) recovers the downwash expression given by Cones.!7 
The general lift circulation distribution is expressed as 
T=P,+T. (66) 


where I, and [, are the rigid and aeroelastic components of the lift circulation distribution which are given 
by 


1 
lr, = 9 Voth, (67) 


(68) 


Be Vast ( 1 — tana® tanaW? °) 


1+tan?a 2 


Thus, the induced angle of attack is not only a function of the lift distribution but also is a function of 
the axial displacement, bending deflection, and torsional twist. 


Consider a special case of the ideal elliptical lift distribution for a rigid wing T, = T94/1—- (2307) 


where « = +U is the coordinate of the elastic axis for the deformed wing and U; = U (8). We use the 
coordinate of the deformed wing to specify lift distribution to ensure that the solution is physical. Otherwise, 
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the lift circulation at 2 = +3 does not exist since x € [-8 — Ui, 8 + U;). Then, the downwash expression 
becomes 


- a (zo - 2’) (1 + Uz.) + (Wo — Oo — W + €0) W. 0 dang 
w; (20) = — 
Eb, am (a — 2)? + (Wo — 8) — W + 0)] 59” 


AT ox. aCe 2 


b2 1-( 22" iy dar’ 


b+2U;, 


For a rigid wing, the ideal elliptical lift distribution produces a uniform induced downwash and the 
minimum induced drag. The effect of the wing bending is to alter this uniform downwash for the ideal 
elliptical lift distribution. As a result, the ideal elliptical lift distribution design for a rigid wing will not 
yield the minimum induced drag due to the nonlinear large bending deflection. Wing twist is normally used 
to change the lift distribution design for a flexible wing. 

As an example, consider a constant-section high aspect ratio wing with 1-ft chord and 30-ft span. The 
cross section is a NACA 0012 with a wall thickness of 0.25 inches and a modulus of elasticity of 10 x 10° 
psi for aluminum. The cross sectional area and the flapwise area moment of inertia are computed to be 


A = 2%kacty and I = §kyctty ) where k4 = 1.02084 and ky = 1.03716. Without considering the 


nonlinear large deflection effect, the wing is designed for an elliptical lift distribution | = 1004/1 — (2) lb/ft 


where | = poo VI for poo = 2.37756 x 10-3 slug/ft® and V,, = = 200 ft/sec corresponding to M,, = 0.2687. 


The lift curve slope is estimated from the formula Cy, = ee where AR = 30 is the aspect ratio. 
co! TA 


This yields a = 6.89°. Figure 6 shows the linear and nonlinear deflections of the wing computed by applying 
the Galerkin method with 10 mode shapes to Eq. (60). Both the wing tip bending deflections as predicted 
by the linear and nonlinear bending theories are in close agreement The nonlinear bending theory yields a 
wing tip deflection of 3.0807 ft or 20.54% of the wing semi-span. In comparison, the linear bending theory 
predicts a wing tip deflection of 3.2137 ft or 21.42% of the wing semi-span. Thus, the nonlinear bending 
theory predicts a smaller bending deflection than the linear bending theory. This difference is due primarily 
to the aeroelastic contribution of the last term in the left hand side of Eq. (60). 


—— Linear Deflection 
—— Nonlinear Deflection 


0 5 10 15 
x+U, ft 


Figure 6. Linear and Nonlinear Bending Deflections 


13 of 26 


American Institute of Aeronautics and Astronautics 


The significant difference between the linear and nonlinear bending theories is the axial displacement. 
For the linear bending theory, the axial displacement is exactly zero since there is no applied axial force. On 
the other hand, the nonlinear bending theory predicts an axial displacement of -0.3584 ft or -2.39% of the 
wing semi-span. As a result, the linear bending theory actually predicts an increase in the wing length by 
2.55% of the wing semi-span, whereas the nonlinear bending theory preserves the length of the wing. This 
length-preserving property of the nonlinear bending theory in effect causes a reduction in the wing aspect 
ratio, which leads to an increase in the induced drag and decrease in lift as compared to the linear bending 
theory. 

The length-preserving property of the nonlinear bending theory can be shown by examining Eq. (56) 
and ignoring the terms ©? and W2,. This leads to 


Uy = \/1— (sin7! W,)? — 1 = -; (sin~! Wy)? —5W2 (70) 


Thus, it can be seen that U, is highly dependent on W,,. As the bending deflection increases, the axial 
displacement decreases. 
The change in length due to the nonlinear bending theory is computed as 


d wa 
pV (1+ Us) + WE = f+ =e (71) 


which shows the length-preserving property of the nonlinear bending theory. 
For the linear bending theory, the length is not preserved. Therefore, the increase in length is computed 


L L 
aL= f (V1+W2 - 1) dic © > if W2dx (72) 
0 0 


This implies that the linear bending theory yields a non-physical solution. 
Also ignoring the term W2, in Eq. (60), the nonlinear static aeroelastic equation for bending can be 
2 
approximated by the Taylor series of sin! W, + W, (1 + me). Then, the resulting equation is expressed 
as 


as 


0? Ww? tanaW? 
Woe ta vo = CL, Goo — m 73 
a | oe ( "% ) | tend “O(i Bian’ a) Pe vo) 


Recognizing that W4 ~ 0, then Eq. (73) becomes 


Oo" tanaW? 


as (E1ly,Wee) + CS ‘50 C — 74 
( yy ) CLA “5 (1 + tan” a) CL, doo€ — Mg ( ) 


Ox? 
Thus, the difference in the bending deflection between the linear and nonlinear bending theories is due 
to mostly the aeroelastic term in the left hand side which contributes positively to the wing stiffness. For 
sweptback wings, truss-braced wings,!® and rotary wings, the effect of the axial tension P, can be a dominant 
factor. Then, the significance of the nonlinear large deflection bending can be much more pronounced. 
Figure 7 is the effective nonlinear elastic angle of attack due to the term W2 in Eq. (54). This shows 
that the nonlinear large deflection bending causes a negative angle of attack of about —0.24°. This results in 
a non-elliptical lift distribution as shown in Fig. 8 which is compared to the ideal elliptical lift distribution 
for the rigid wing and for the flexible wing using the linear bending theory. To illustrate the effect of length- 
preserving property of the nonlinear bending theory, the lift distribution for the nonlinear large deflection 
bending is plotted against the displaced coordinate. To restore the ideal elliptical lift distribution, the wing 
would have to be re-twisted along the wing span to zero out the effect of the elastic angle of attack. Figure 
8 shows the elliptical lift distribution plotted in the displaced coordinate. 
Figure 9 is the plot of the induced angle of attack a; = 4 due to the induced downwash over the wing 


Voc 
as computed by Eq. (69) which, for the constant chord and negligible torsional twist, is expressed as 


(zo) ie (x - 2’) (1+ Uz,0) + (Wo — Oo — W + €0) We. AD ox" ellie 
Pita t 7 v 
abou; 4n |(2 x')” + (Wo — eo — W 4 0)"| bef 1 — (230 ) . 
+2U;¢ 
(75) 
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which takes advantage of the length-preserving property of the nonlinear bending theory. 


X, ft 


Figure 7. Aeroelastic Angle of Attack due to Nonlinear Bending Deflection 
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— Elliptical Distribution for Nonlinear Deflection with Re-Twist 
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2(x+Uyi(b+2U) 


Figure 8. Lift Distributions for Bending Deflections 


Note that for the linear bending theory, [. = 0 and U = 0. So, the induced downwash is computed from 
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the following equation: 


ilo 


wise (x9) = / (Xo x) t (Wo eOo W + eQ) Wz,o 1 Tox a 
-3 (xo — 2)? + (Wo —eOp —-W + 0)"| V2 +W?2, be/1 _ (4) 


(76) 


The rigid wing has a constant induced angle of attack for the ideal elliptical lift distribution whereas the 
bending deflection causes the induced angle of attack to become non-uniform with the maximum value at 
the mid-wing span. The linear bending theory results in a lower induced angle of attack than the nonlinear 
deflection. This is due to the difference in the wing span where for the nonlinear deflection there is a reduction 
in the wing span which causes an increase in the induced downwash, hence induced angle of attack. The 
induced angle of attack at the center of the wing is the greatest with the nonlinear deflection. The results 
illustrate the significance of the effect of the nonlinear bending theory on the induced drag of a wing structure. 
Also shown in Fig. 9 is the induced angle of attack for the nonlinear deflection when the wing is re-twisted 
to remove the effect of the nonlinear elastic angle of attack. This has an effect of lowering the maximum 
induced angle of attack for the nonlinear deflection. 


0.51 


— Linear Deflection 
—— Nonlinear Deflection 
—— Nonlinear Deflection w/ Re—Twist 
— Rigid 


-15 -10 -5 


0 
x+U, ft 


Figure 9. Induced Angles of Attack due to Bending Deflections 


The lift coefficient is evaluated by taking into account the wing bending deflection slope and the shortening 
of the wing span due to the axial displacement according to 


b b 
1 2 1 atl , 
C..= as Poo VooI cos (sin~* W,,) (1+ Uz) dz = aS PooVool V1 — W2dx (77) 
fore) -3 oe) —3-U, 


The induced drag coefficient is evaluated from the following equation: 


1 $4U; 
a. 


a 


dooS —$-U; 


Table 1 shows the lift and induced drag coefficients and the span efficiency factors € based on the aspect 
ratio of the rigid wing for the linear deflection, nonlinear deflection, and rigid wing. The lift and induced drag 
coefficients are evaluated with the reduced wing span by 2U; for the nonlinear deflection. It can be seen that 
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the effect of the nonlinear deflection is to reduce both the lift coefficient and the span efficiency. The linear 
bending theory predicts only a 2.34% lift reduction and 1.19% decrease in the span efficiency relative to the 
rigid wing. In comparison to the rigid wing, the nonlinear bending theory predicts a 6.69% lift reduction 
and 9.02% decrease in the span efficiency. Even with re-twisting, the lift reduction only decreases to 4.88%. 
Thus, the linear bending theory tends to yield more optimistic aerodynamic performance characteristics for 
a high aspect ratio flexible wing. 


Cr Cp; € 
Linear Deflection 0.7166 | 0.005514 | 0.9881 
Nonlinear Deflection 0.6847 | 0.005468 | 0.9098 
Nonlinear Deflection with Re-Twist | 0.6980 | 0.005515 | 0.9375 
Rigid Wing 0.7338 | 0.005716 1 


Table 1. Lift and Drag Coefficients and Span Efficiency Factors 


V. Weak-Form Numerical Solution 


The nonlinear partial differential equations are usually referred to as a strong form. A numerical so- 
lution can be performed by converting the strong-form partial differential equations into weak-form in- 
tegral equations by letting U (#,t) = Nu (x)u(t), V(a,t) = Ny (xz) v(t), W(2,t) = Nw («) w(t), and 
O (x,t) = No (x) @(t). For the Galerkin method, Ny, (x), Nu (x), Nw (x), andNg (x) are the vectors of the 
mode shapes over the domain of the solution x € [0,Z], whereas for the finite element method they are 
the vectors of interpolation functions over the domain of the solution for an element x € [0,1]. The vec- 
tors u(t), v(t), w(t), and @(t) represent the generalized coordinates in the Galerkin method, or the nodal 
displacements including the bending slopes in the finite element method. 


2 
The term sin~! W, can be approximated by the first two terms of the Taylor series expansion W, (1 + Le ) s 


Pre-multiplying Eqs. (30) to (33) by N/ (x), N} (x), NJ, (x), and Nj (a), respectively, and integrating over 
the domain of the solution yields 


; 1 Life 
= [Ncza (1 be 50) Ni udr — 5 ine 


W? : 7 Were at 
EAW, (1+ - ) Ni, 4 El tie wate 


1 q , 
25 . N,, ElzaQ2N,0dx = — / N) frdx (79) 


= i N," (GJ + Prr?) N,Odx = — ‘i Nj medx (80) 
uw 1 a Uz ” , , om 


wt wt / / if aw 2 t 
i (NUTEL.N, 4 N,"P.N,,) vdx — 5 | Nv EL.®: (1 fe “) N.,wdx 


- Om: 


1 uw / 
- 5 [ No" Els sin”' W,N,0dx =f Ny (4, 5 ) ae (82) 
x 


The problem is cast as 
K (U,W,0)x =F (83) 


om 
where x = u’_ v' wi @! , K is the nonlinear stiffness matrix, and F is the force vector. The 


nonlinear global stiffness matrix and force vector are assembled from the following elemental stiffness matrix 
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and force vector: 


kyu 0 kuw ku 
k= 0 kyy kuw ku (84) 
0 O kyy 0 
0 0 0 koe 
Ni fc 
Ny (fy — one ) 
f= me ye pa 
i i (j, + 8m) | (85) 
Ny Mx 
where : 
oe / NBA (1 £ 50») Ni de (36) 
k =5/NI raw, (14 82) Ni, 4 21, =n" ae (87) 
uw — 2 uU x 6 wi! wyy = Ww2 w 
1 a fi 
ku = 5 [| Ny! ElrxO2Node (88) 
i= if (Ni EI,,N,dx + N," PsN.) da (89) 
1 u” Ww? , 
how =—5 f NoTBlnOs (14+ 42) Nude (90) 
1 w" / 
kvo = 5 [NOEL sin”! W,Nogdax (91) 
uw 1 £ ” , , 
Kow = iF A An a5 N,, PrN, dx (92) 
1_W2 
ine ‘i Nj" (GJ + Per?) Node (93) 


Note that the flapwise bending stiffness includes both the geometric nonlinear large deflection effect due 


to the term a and the nonlinear axial tension stiffening effect due to P, which also exists in the 


torsional stiffness. The nonlinear large bending deflection effect thus increases the flapwise bending stiffness 
which effectively reduces the bending deflection. 

The displacement matrix equation is nonlinear and can be solved for the static aeroelastic deflections by 
any nonlinear root search methods. One iterative method is proposed as follows: 


e For i = 0, initialize U; (x) = 0, V; (x) = 0, W; (x) = 0, and ©; (x) = 0. Evaluate K;. Then, compute 
x; =K;'F (94) 
e For i =1,2,...,n, evaluate U; (x), V; (x), Wi (a), and ©; (x). Update K;. Then, iterate 
x; = mx;_1 + (1—m)K;'F (95) 


with 0 < m < 1 until x; converges. Note that this algorithm will ensure numerical stability and 
converge to the ’exact’ solution of Eq. (94) in the limit as n + oo. Equation (95) essentially imposes 
the stability condition on the unit circle. Otherwise, if Eq. (94) is used instead which corresponds to 
m = 0, numerical stability issues sometimes may be encountered. A suitable value for m is usually in 
the range between ; and 5. 
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VI. Computational Results 


A coupled aerodynamic-nonlinear finite element model is constructed for the CRM sub-scale model. 
The flow condition is at Mach 0.1162 and a dynamic pressure of 20 psf. The model is developed within a 
multidisciplinary coupled aerodynamic-nonlinear finite element analysis framework as shown in Fig. 10. The 
framework consists of four computational modules: 1) NASA vortex lattice code VORLAX?® developed by 
Louis Miranda et al. with a pre-processor VORVIEW for providing a graphical user interface (GUI) for model 
geometry handling, 2) 2D Euler CFD code MSES with an integral boundary layer method developed by Mark 
Drela,?° 3) an in-house computational geometry code that manipulates the model geometry deformation, 
and 4) an in-house nonlinear 3D beam finite element analysis (FEA) code with structural dynamics and 
aeroelasticity capabilities. The framework can compute the solution for a specified lift coefficient input 
or angle of attack input. For each solution for lift matching, VORLAX iteratively computes the angle of 
attack to match the specified lift coefficient. Then, the solution is corrected for the viscous and, if necessary, 
transonic flow effect with MSES which computes the lift and drag coefficients for each wing station in the 
VORLAX model. This transonic and viscous correction method requires an iteration loop between MSES 
and VORLAX. Once the corrections are made, VORLAX determines the new angle of attack iteratively to 
match the specified lift coefficient. At this point, the aerodynamic forces and moments are computed and 
then passed to the FEA code. The FEA code computes the wing deflections and then passes this deflection 
information to the computational geometry module which computes the geometry of the model deformation. 
This updated geometry is then passed to VORVIEW for geometry processing. The solution process is then 
repeated until the angle of attack converges to within a specified tolerance. 


Lift 

Coefficient: 

Cc, Geometry Vorview/ VORLAX 
VORLAX 


Lift matching loop Lift matching loop 


Estimate a 


Estimate a 


Transonic and viscous correction loop 


Aeroelastic loop 


Figure 10. Coupled Aerodynamic-Nonlinear Finite Element Analysis Framework 


Figure 11 shows the convergence of the iterative nonlinear solution according to Eq. (95) for three values 
of m = 0, 5, 5. The wing tip deflection W; at the angle of attack of 7.616° is plotted as a function of 
the iteration number. The results show that the solution converges the fastest with m = 0 after only two 
iterations. As m increases, the convergence is slower. It is noted that in other situations where the nonlinear 
effects are much more dominant, experiences support the use of a positive value for m between 3 and 5. 

Figure 12 shows the linear and nonlinear bending deflections of the CRM sub-scale model along the 
aircraft pitch axis y computed for the angle of attack of 9.538°. The linear deflection at the wing tip is 10.17 
inches versus the nonlinear deflection of 10.11 inches. The wing length measured along the elastic axis is 
99.69 inches. So the wing tip deflection is about 10% of the wing length. The nonlinear deflection is only 
slightly smaller than but essentially the same as the linear deflection. 

Figure 13 shows the linear and nonlinear axial displacements of the CRM sub-scale model along the 
aircraft pitch axis y. The linear theory predicts an axial displacement of 0.01456 inches at the wing tip 
whereas the nonlinear theory predicts an axial displacement of -0.7124 inches. The axial extension with 
the linear theory is due to the drag force component along the swept elastic axis. As expected, the axial 
displacement computed by the linear theory is non-physical since it predicts a larger elongation of the wing 
than is possible. The linear theory predicts an elongation of 0.8396 inches as compared to the nonlinear 
theory which predicts an elongation of only 0.1136 inches due to the drag component along the elastic axis. 

Figure 14 shows the linear and nonlinear chordwise bending deflections of the CRM sub-scale model along 


the aircraft pitch axis y. As expected, the chordwise bending deflections are very small due to the much 


19 of 26 


American Institute of Aeronautics and Astronautics 


larger chordwise bending stiffness and the small drag force component acting normal to the elastic axis. 
Nonetheless, the difference between the linear and nonlinear bending theories is revealing. The nonlinear 
bending theory predicts a larger chordwise bending deflection but in the opposite direction than the linear 
bending theory. This is due to the contribution of the nonlinear term sin~' W,,0, in Eq. (33) which causes 
a negative contribution to the chordwise bending deflection. 
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Figure 11. CRM Sub-Scale Model Solution Convergence at a = 7.616° 
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Figure 12. CRM Sub-Scale Model Flapwise Bending Deflection at a = 9.538° 
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Figure 13. CRM Sub-Scale Model Axial Displacement at a = 9.538° 
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Figure 14. CRM Sub-Scale Model Chordwise Bending Deflection at a = 9.538° 


Figure 15 shows the linear and nonlinear torsional twists of the CRM sub-scale model along the aircraft 
pitch axis y. Both the linear and nonlinear torsional twists are virtually the same. The nonlinear theory 
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predicts a smaller torsional twist than the linear theory due to the presence of the small axial tension caused 
by the drag component along the elastic axis which acts to increase the torsional stiffness. This is observed 
in the results which show that the nonlinear torsional twist at the wing tip of —0.5502° is slightly smaller 
than the linear torsional twist of —0.5510°. 
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Figure 15. CRM Sub-Scale Model Torsional Twist at a = 9.538° 
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Figure 16. CRM Sub-Scale Model Flapwise Bending Deflection vs. Bending Stiffness 
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Figure 16 shows the effect of the flapwise bending stiffness on the linear and nonlinear bending deflections 
at the wing tip. The nonlinear bending deflection is generally smaller than the linear bending deflection, 
but as the stiffness is reduced by 50%, the trend is reverse. As the stiffness increases, the nonlinear and 
linear bending deflections appear to converge. This is expected because as the stiffness increases the bending 
deflection decreases until it is small enough that the linear bending theory holds. 
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Figure 17. CRM Sub-Scale Drag Model Lift Coefficient vs. Angle of Attack at Mj = 0.1162 
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Figure 18. CRM Sub-Scale Drag Model Lift Coefficient at a = 7.616° and Mx = 0.1162 vs. Bending Stiffness 
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Figure 20. CRM Sub-Scale Model Drag Polar 


Figure 17 shows the lift curves of the CRM sub-scale model. The rigid wing produces more lift at a given 
angle of attack than the flexible wing as expected. Both the linear and nonlinear bending deflections give 
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almost the same lift curves. On a closer look, the lift curve for the linear bending deflection is slightly higher 
than that for the nonlinear bending deflection. This becomes more apparent in Fig. 18. As the wing becomes 
stiffer, the difference in the lift coefficient between the nonlinear and linear bending deflection decreases. 
Figure 19 shows the drag polar of the CRM sub-scale model. For the wind tunnel test, a jig shape twist 
optimization has been conducted to determine the optimal jig-shape twist that minimizes the induced drag 
for the CRM sub-scale model in its deflected state in the wind tunnel. As a result, at a fixed value of the lift 
coefficient, the rigid wing produces a higher drag than the flexible wing. The linear and nonlinear bending 
deflections both produce nearly identical drag coefficients. On a closer look, the drag coefficient for the 
nonlinear bending deflection appears to be slightly smaller than that for the linear bending deflection but 
not by much. At the angle of attack of 8°, the drag coefficients for the linear and nonlinear deflections are 
0.020560 and 0.020557, respectively. Figure 20 shows that the drag coefficient for the nonlinear bending 
deflection is slightly lower than that for the linear bending deflection as the wing bending stiffness varies. 
Figure 21 is the color map of the surface pressure distribution on the CRM sub-scale model. 
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Figure 21. CRM Sub-Scale Model Surface Pressure Map at a = 7.616° and Mx = 0.1162 


VII. Conclusions 


This paper presents a nonlinear bending theory for large deflection bending. The nonlinear bending 
theory properly accounts for the length preservation when the structure is not subject to axial loading. 
While the linear bending theory predicts a similar bending deflection characteristic, it does not account 
correctly for the axial displacement. As a result, the linear bending theory is not length-preserving. At 
smaller bending deflections, the length-preserving property may not be significant, but when the bending 
deflection is sufficiently large, the nonlinear bending theory predicts a reduction in the wing span whereas 
the linear theory would fail to predict this effect. As a result, this could cause a significant effect on the wing 
aerodynamic performance. A modified lifting line theory is developed to take into account the full deflection 
states of the wing deformation. An analysis is conducted via an example of a high aspect ratio unswept wing 
to illustrate the effect of nonlinear large deflection bending on aerodynamics based on the lifting line theory. 
The results show that the nonlinear bending deflection causes a reduction in lift and span efficiency factor. 
The induced downwash for the ideal elliptical lift distribution cannot remain uniform across the wing span 
due to the bending deflection. 
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A coupled aerodynamic-nonlinear finite element model is developed for a Common Research Model 
(CRM) sub-scale wind tunnel model that will be tested at the University of Washington Aeronautical Lab- 
oratory to demonstrate an active real-time drag optimization control strategy. The computational results 
show that for this particular wind tunnel model there is a slight reduction in lift with the nonlinear bending 
theory, but otherwise the nonlinear bending theory does not have a dominant effect on the aerodynamic 
performance as compared to the linear bending theory. This study suggests that the linear bending theory 
is a good approach for analyzing this wind tunnel model. 
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